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Sobre este documento 


El presente documento que se adjunta al visor web de susceptibilidad a incendios forestales del 
cantón Loja pretende ser un complemento técnico para su comprensión y utilización. En él tendrá 
una vista rápida del proceso de recolección de información, entrenamiento y validación de 
modelos de aprendizaje automático, así como, la automatización en la plataforma Google Earth 
Engine con una explicación del proceso para lanzar una Earth Engine App y un glosario de 
términos utilizados. 
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IDENTIFICACIÓN DE VARIABLES EXPLICATIVAS 


La generación de la base de datos para el cantón Loja tuvo la finalidad de recolectar información 
de variables relacionadas con el inicio de incendios forestales (momento de ignición), a partir de 
una revisión de literatura científica. Esta base de datos permitirá generar modelos de 
susceptibilidad de incendios forestales para identificar zonas con una probabilidad de ocurrencia 
mayor y que, por tanto, necesitan una atención prioritaria en la prevención de estos desastres. 


Este apartado del documento surge en colaboración con tres estudiantes del componente 
práctico Gestión Productiva 3.1. de la titulación de Gestión Ambiental de la Universidad Técnica 
Particular de Loja, y pretende dar una breve explicación del proceso de obtención de las 23 
variables que conforman la base de datos. En función de ello, esta sección se subdivide en: 
variables topológicas, variables de interacción antrópica y variables del estado de la vegetación. 


VARIABLES TOPOGRÁFICAS 


Variable Descripción Proceso de obtención | Metadatos de la capa 
Se importa un archivo .kml con los 
límites del cantón de interés y se 
realiza la búsqueda de Modelos de 
Elevación Digital (DEM) disponibles en А е 
la web de ALOS PALSAR para Іа il ias bl 
zona; Resolución espacial: 12.5 x 
Una vez descargados los archivos en 12:3 metros 
high resolution se utiliza la А 
1 Elevación Altura sobre el herramienta de QGIS para combinar о valores: 
nivel del mar. las teselas (Ráster, Miscelánea, 
Combinar), posteriormente se recorta Unidades mñn 
a los límites de nuestra zona de ages: тзг; 
interés (Ráster, Extracción, Cortar 
ráster por capa de máscara) 
SRC: EPSG:32717 
AI DEM obtenido en la web de ALOS a ЕТЕ 
РАІ ЅАВ, se le aplica un filtro de Resölución espacial: 12.5:х 
Grado de - _ 12.5 metros 
inclinació suavizado en QGIS herramienta 
inclinación del SAGA (Simple Filter, S th3 
eee puede (Simple Fiton Smoot g), Rango de valores: 
. medirse en - 
2 Pendiente grados o Después de aplicar el filtro de 0-375 
¡ suavizado, se utiliza la herramienta a : 
а GRASS (r.slope.aspect, percent) Unidades: porcentaje (%) 
Dirección de 
pendiente SRC: EPSG:32717 
descendente Al DEM obtenido en la web de ALOS Е А 
que presenta PALSAR, se le aplica un filtro de ree espacial: 12.5 x 
{ кле nc suavizado еп QGIS herramienta metros 
omando como SAGA (Simple Filter, Smooth 3). 
| Е referencia los (Simp ) Rango de valores: 
Е оеп Pard А Después de aplicar el filtro de а 
suavizado, se utiliza la herramienta А . o 
(Norte, Sur, GRASS (r.slope.aspect) di 
Este, Oeste) 


VARIABLES DE INTERACCIÓN ANTRÓPICA 


La generación de las capas de accesibilidad surge a partir de una capa de superficie de 
fricción,que representa el costo añadido que implica la movilización desde cada punto del 
territorio hasta el punto(s) objetivo(s), la siguiente gráfica pretende ilustrar el proceso seguido 


para su cálculo. 


Centro Económico (4) 


Obtenido en Mapa Interactivo Ambiental Generado en SAGA, apartir del DEM Obtenido en IGM (hrtps://sná.gob.ec/ Generado en QGIS, apartir del DEM Áreas Protegidas (5 
(btrp:/fide.ambiente.gob.ec/ ALOS PALSAR coberturas), actualizado con digitalización ALOS PALSAR ES 
) (hups://search,asf,alaska.edu/5/) manual (rrps://search,asf.alaska.eduls/) 
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Metodología tomada de: 


Reyes, F. (2017). Mapa de Accesibilidad Territorial. En Manual de Prácticas de Ordenación Territorial (págs. 29- 
41).Loja: UTPL. 


Variable Descripción Proceso de obtención | Metadatos de la capa 
Accesibilidad rio: AR ¡teams la herramienta аа 
сөбө cualquier punto del | GRASS r.cost.full, se configura esolución espacial: 10 
económico а а е! сото рипїо(в) objetivo(s) la oda 1 Е 
4 сепіго есопбтісо сара Че centro:económico 0 aA уи 
` | descargada en el portal del SNI. | Unidades: minutos 
Utilizando la herramienta 
GRASS r.cost.full, se configura 
como punto(s) objetivo(s) la 
capa de áreas protegidas 
descargada en el Mapa | 
Accesibilidad Facilidad de кезиин. A a 
a áreas movimiento desde А f 10 х 10metros 
5 protegidas cualquier punto del | Las categorías consideradas de | Rango de valores: 
cantón hasta áreas | aca Mona de ioes Protegias | ОА дев: 
pr otegidas. Bosque y Vegetación Protectora Unidades: minutos 
Area Bajo Conservación 
Patrimonio Forestal del Estado 
Areas de Protección Hídrica 


SRC: EPSG:32717 


Facilidad de Resolución espacial: 10 
Accesibilidad a movimiento desde Utilizando la herramienta x 10metros 
6 | zonas pobladas | cualquier punto del | GRASS r.cost.full, se configura a Че valores: 
cantón hasta como punto(s) objetivo(s) la + є 
poblados. capa Кү А к; Unidades: minutos 
en el portal del SNI. 
SRC: EPSG:32717 
Facilidad de Utilizando la herramienta Resolución espacial: 10 
Accesibilidad a movimiento desde GRASS r.cost.full, se configura | х 10metros 
7 vías cualquier punto del como punto(s) objetivo(s) la pango de valores: 
cantón hasta vías. capa de vías descargada en el Ж Л 
йн, del SNI y a ABLA RS minutos 
manualmente. 
Facilidad de Utilizando la herramienta SRC: EPSG:32717 
movimiento desde GRASS r.cost.full, se configura | Resolución espacial: 10 
Accesibilidad a cualquier punto del como punto(s) objetivo(s) la x 10metros 
8 zonas cantón hasta zonas capa de cobertura vegetal Rango de valores: 
antrópicas pobladas o zonas descargada en el Mapa 0 - 262 f 
agropecuarias. Interactivo Ambiental. Unidades: minutos 
SRC: EPSG:32717 
Distancia Utilizando la herramienta Resolución espacial: 10 
desde el GRASS r.grow.distance, se x 10metros 
9 centro Distancia desde el configura como punto(s) Rango de valores: 
económico centro económico. objetivo(s) la capa de centro 0- 4173 
económico descargada en el Unidades: metros 
portal del SNI. 
SRC: EPSG:32717 
Distancia . | Utilizando la herramienta Resolución espacial: 10 
desde ríos Distancia desde GRASS r.grow.distance, se | х 10теігоѕ 
10 ríos. configura como punto(s) Rango de valores: 
objetivo(s) la capa ríos generada | 0- 4173 
a partir del DEM. Unidades: metros 
. | Utilizando la herramienta 
Distancia Distancia desde GRASS r.grow.distance, se SRC: EPSG:32717 
desde vías vías. configura como punto(s) Resolución espacial: 10 
objetivo(s) la capa vías х 10metros 
11 descargada en el portal del SNI | Rango de valores: 
y actualizada manualmente. 0— 14253 
Unidades: metros 
Utilizando la herramienta 
Distancia GRASS r.grow.distance, se SRC: EPSG:3271 7 
12 desde áreas Distancia desde configura como punto(s) Resolución espacial: 10 
protegidas áreas protegidas objetivo(s) la capa áreas х 10metros | 
protegidas descargada en el апе valores: 
Mapa Interactivo Ambiental. Unidades: metros 
Las categorías consideradas de 
áreas protegidas son: 
Sistema Nacional de Áreas Protegidas 
Bosque y Vegetación Protectora 
Area Bajo Conservación 
Patrimonio Forestal del EstadoÁreas 
de Protección Hídrica 
SRC: EPSG:32717 
Utilizando la herramienta Resolución espacial: 10 
Distancia a Distancia desde GRASS r.grow.distance, se x 10metros 
13 zonas zonas pobladas o configura como punto(s) Rango de valores: 
antrópicas zonas objetivo(s) la capa zonas -661 
agropecuarias. pobladas o zonas Unidades: metros 


agropecuarias. 


VARIABLES DEL ESTADO DE VEGETACIÓN 


i F Proceso de 
Variable Descripción obtención Metadatos de la capa 
Atmospherically Resistant 
Vegetation Index 
NDVI corregido para efectos de SRC: EPSG:32717 
ispersión atmosféri 
14 ARVI dispersión atmosférica Resolución espacial: 10 x 10 
arvi „ МЕС (2 + Red) + Blue) metros 
NIR — (2 * Red) + Blue) Rango de valores: -1 a 1 
NIR (Near-Infrared Region) > B8 
Red (Rojo, espectro visible) > B11 
Blue (Azul, espectro visible) > B2 
Bare Soil Index 
Identificar zonas sin cobertura vegetal 
poj = Bed +SWIR)=(NIR + Blue) SRC: EPSG:32717 
(Red + SWIR) + (VIR + Blue) 2. 
Los Índices Resolución espacial: 10 x 10 
15 BSI Red (Rojo, espectro visible) > B11 | espectrales o índices | metros 
SWIR (Shortwave Infrared) > B4 А 
NIR (Near-Infrared Region) > B8 de vegetación son R d | :-1а1 
Blue (Azul, espectro visible) > B2 descargados ango:06 valores: =la 
mediante el siguiente 
Green Normalized Difference | script de Google Earth 
Vegetation Index Engine, se иза para 
ello imágenes Sentinel 
A , 2 2А. А А 
NDVI más sensible a la variación del SRC: EPSG:32717 
16 GNDVI contenido de clorofila y nitrógeno Resolución espacial: 10 x 10 
сару 2 MR- SWIR тое 
NIR + SWIR Rango de valores: -1 a 1 
SWIR (Shortwave Infrared) > ВЗ 
NIR (Near-Infrared Region) > B8 
Moisture Stress Index 
SRC: EPSG:32717 
i7 MSI Determinar el eden щш de la Resolución espacial: 10 х 10 
vegetació metros 
ms] = MidiR Rango de valores: 0 a 3+ 
NIR 
MidIR (Shortwave Infrared) > B11 
NIR (Near-Infrared Region) > B8 
Normalized Burn Ratio 
Estimar la cantidad y desarrollo de SRC: EPSG:32717 
areas:incendiadas Resolución espacial: 10 x 10 
18 NBR _ NIR — SWIR 110x 
NBR = MIR + SWIR metros 
SWIR (Shortwave Infrared) > B12 Rango de valores: -1 a 1 
NIR (Near-Infrared Region) > B8 
Normalized Difference Moisture 
Index SRC: EPSG:32717 
19 NDMI Determinar и agua de la Resolución espacial: 10 x 10 
NIR — SWIR metros 
NDMI = JIR + SWIR R de valores: -1 a 1 
SWIR (Shortwave Infrared) > B11 апдо бө магогез ка 
NIR (Near-Infrared Region) > B8 
Normalized Difference Water 
Index SRC: EPSG:32717 
20 NDWI e zonas de Resolución espacial: 10 x 10 
NDW[= NIR — Green metros 
NIR + Green Rango de valores: -1 a 1 
SWIR (Shortwave Infrared) > B3 
NIR (Near-Infrared Region) > B8 
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МОМ! 


Normalized Difference 
Vegetation Index 
Facilidad con la que el fuego se 
propaga en la vegetación 

NIR — Red 


NDVI = те Red 


Red (Rojo, espectro visible) > B4 
NIR (Near-Infrared Region) > B8 
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SAVI 


Soil Adjusted Vegetation Index 
Analizar la vegetación en etapas de 
crecimiento inicial o dispersa 


NIR — R 
* (1+ ZL) 


SAVI = NIR+R 4L 


NIR (Near-Infrared Region) > B8 
Red (Rojo, espectro visible) > B4 
Factor de correción de la luminosidad del suelo > L 


23 


Moisture 
Index 


Normalized Difference Moisture 
Index 
Estimar el contenido de humedad de 


la vegetación 
NIR — SWIR 


Moisture = WIR 


SWIR (Shortwave Infrared) > B11 
NIR (Near-Infrared Region) > B8A 


SRC: EPSG:32717 


Resolución espacial: 10 x 10 
metros 


Rango de valores: -1 a 1 


SRC: EPSG:32717 


Resolución espacial: 10 x 10 
metros 


Rango de valores: -1 a 1+ 


SRC: EPSG:32717 


Resolución espacial: 10 x 10 
metros 


Rango de valores: -1 а 1+ 


GENERACIÓN DE LA BASE DE DATOS 


Una vez culminada la identificación y espacialización de las variables, se dio paso a la generación 
de la base de datos. En este proceso se utilizó el software libre QGIS v.3.22. con su complemento 
Point Sampling Tool, el que permite capturar valores de capas ráster o vectoriales que les 
corresponden a los puntos de muestreo (puntos de ocurrencia y no ocurrencia). De las 23 variables 
contenidas en la base de datos, 13 tienen un carácter permanente y 10 son actualizables 
(disponibilidad de imágenes Sentinel 2A cada 5 días), por lo que el tratamiento de estos datos es 


distinto. Para ilustrar mejor este proceso se puede consultar el siguiente diagrama de flujo: 


formato ráster (TIFF 


5 capas de distancia 


formato ráster (TIFF) 


| 5 capas de accesibilidad 


3 capas de topografía 


formato ráster (TIFF) 


10 capas de índices de vegetación 


formato ráster (TIFF) 


(actualizable) 


Seleccionar por expresión (fecha) y 
actualizar tabla de datos (solo objetos 


seleccionados) 


Extraer valores con Point Sampling Тоо! 


| 


formato vectorial (.shp 


Base de datos en formato vectorial (.shp) 


ма capa de puntos de muestreo 


Base de datos en 
formato .csv 


ENTRENAMIENTO Y VALIDACIÓN DE LOS MODELOS ОЕ SUSCEPTIBILIDAD A 
INCENDIOS FORESTALES 


La generación del modelo de susceptibilidad siguió el proceso que se explica en la siguiente figura: 


APRENDIZAJE AUTOMÁTICO 


Es una rama de la inteligencia los sistemas pueden aprender de datos, 
ы identificar patrones y tomar decisiones con 
artificial (IA) e 
mínima intervención humana 


| ¿qué se desea hacer? Reunir datos de fuentes Limpiar y dar formato | 
Problema a resolver confiables 
E ЭЛИ 
Entrenar el Validar el modelo 


Proceso de Machine Learning 


modelo 
Técnicas de aprendizaje: н Evaluar los resultados Ingresar nuevos 
-  Supervisado Corer e redor mo con los datos restantes datos en el modelo 
-  Nosupervisado datos disponibles (25%) a confiable 
l › 1.Árboles de decisión. 

2. Clasificación de Naive Bayes. 5.Support Vector Machines (SVM). 

3.Regresión por mínimos cuadrados. 6.Métodos “Ensemble” (Conjuntos de 

4.Regresión Logística. clasificadores). 


Ejemplos de técnicas de aprendizaje automático 


(Combustible `\ 
(Índices de + 
vegetación) 


14. ARVI 


Incendios históricos Variables + Cercanía a actividades antrópicas 
(2018-2020) topográficas 
Puntos de ocurrencia (50 %) 1. pendiente Accesibilidad a Distancia desde 


5 2. orientación j 
Puntos de no ocurrencia (50 %) 3, elevación (minutos) (metros) 


9. centro económico 


$ 10. rios 

Satélites VIIRS 3 áreas protegidas ы А áreas protegidas 
MODIS . " . vias . vias 

Ц MOT. завел, 8. zonas antrópicas 13. zonas antrópicas 


sensor PALSAR 


Datos de Datos de resolución: 12.5 m 
validación entrenamiento 


(25 %) (75%) capas base; 
IGM 


SUIA 
SIGTIERRAS 


RECOPILACIÓN DE DATOS Y PREPROCESAMIENTO 


D z 

ШЕ 

El á 

33 

85 М d 
lapas de 

3 Z susceptiblidad a 

о 3 Incendios 

Y Forestales 

De 

ez 

аш 


С) Fuente de descarga Ф Obtenido en Google Earth Engine (О) Procesos [MI] Algoritmos en acis 


Se realiza una partición de las bases de datos en dos grupos: 75% para entrenar el modelo y 25% 
para validarlo según su desempeño en las métricasde evaluación. Para el cantón Loja se probaron 
tres algoritmos: árboles de decisión logística (LMT), regresión logística (LR) y splines de regresión 
adaptativa multivariadas (MARS). 


Por último, en las etapas 6 y 7 del proceso de aprendizaje automático se evalúan los resultados 
con los datos de validación (25% del conjunto de datos total) y, si las métricas son buenas, se 
procede a la predicción con nuevos datos. 


Independientemente, de la técnica o algoritmo empleado, los modelos fueron evaluados con dos 
herramientas: el Area Bajo la Curva (AUC) y el Coeficiente Kappa de Cohen. 


1.1.1. Агеа Bajo la Curva 
La curva ROC (Receiver Operating Characteristive) permite evaluar la capacidad del 
modelo entrenado para clasificar información. Este mide la relación entre la 
sensibilidado “tasa de verdaderos positivos” y la 1-especificidad o “tasa de falsos 
positivos” (Fórmulal y 2); 

Sensibilidad =——— ; VP = verdaderos positivos FN = falsos negativos 


(VP+FN) (1 


Especificidad = YN VN = verdaderos negativos FP = falsos positivos 


(FP+VN) 


Así, el área que se encuentre bajo esta curva representa el porcentaje de 
clasificaciones correctas. Generalmente, los valores superiores a 0.70 son 
considerados como buenospredictores!. El AUC para los modelos entrenados рага el 
cantón fue superior a 0.85. 


1.1.2. Coeficiente Kappa de Cohen 

El coeficiente Kappa o índice Kappa es un instrumento estadístico que estima la 
concordancia entre dos observadores o herramientas de medición para un mismo 
fenómeno, con el propósito de ajustar el efecto del azar en la proporción de la 
concordancia esperada, tal como establece la fórmula 3 ?. 


Kappa = Pate; 


) 


1-Р, 


Р, = proporción de concordancia observada 


Pe = proporción de concordancia esperada por el azar 


Para evaluar la concordancia de variables categóricas, especialmente aquellas con 
salida dicotómica o binaria (como en el presente estudio), Landis y Koch proponen la 
interpretación cualitativa más conocida. En ella, los valores menores a 0.40 
representanuna concordancia pobre o débil; entre 0.41 a 0.6, moderada; entre 0.61 a 
0.8, concordancia buena o significante y aquellos valores mayores a 0.80 tienen un 
nivel deconformidad casi perfecto. 

En los modelos probados para el cantón, el valor del índice Kappa rodea el 0.82. 


1 Cerda, Jaime, 8 Cifuentes, Lorena. (2012). Uso de curvas ROC en investigación clínica: Aspectos teórico-prácticos. 
Revista chilena de infectología, 29(2), 138-141. https://dx.doi.org/10.4067/50716-10182012000200003 

2 Landa, F., Cruz, N., Coria, G., Rebolledo, G., & Hernández, S. (2014). El Índice Kappa de Cohen como Herramienta para 
la Reducción de Ruido de Información Cerebral en la ldentificación de Niveles de Valencia. Puebla. 


AUTOMATIZACIÓN DE LOS MODELOS EN LA PLATAFORMA DE GOOGLE 
EARTH ENGINE 


El diseño y funcionamiento del visor web de susceptibilidad a incendios forestales para el cantón 
Loja está basado en el lenguaje de programación Javascript, es decir, se configura todo desde 
el editor de código de la plataforma de Google Earth Engine. 


A continuación, encontrará una explicación detallada de la implementación del modelo MARS 
para el cantón Loja, si desea tener un mejor entendimiento del fundamento estadístico del 
algoritmo, recomendamos las siguientes lecturas: 


e MARS: Multivariate Adaptive Regression Splines — How to Improve on Linear Regression? | by 
Saul Dobilas | Towards Data Science 


e Chapter 7 Multivariate Adaptive Regression Splines | Hands-On Machine Learning with В 
(bradleyboehmke.github.io) 


Explicación del script para el modelo de susceptibilidad a incendios forestales, 
cantón Loja 


httos://bit. ly/3TaazQ0O 


Líneas 5-15: declaración de variables permanentes (ej: accesibilidad al centro económico, 
distancia desde ríos, distancia desde vías, distancia desde zonas antrópicas, elevación, límite 
cantonal, límites parroquiales) 


B7 = ee.Image("users/lojanjulia/comparacionmodelosIF/Elevacion_div_1000"), 
B6 = ee.Image("users/lojanjulia/comparacionmodelosIF/dist_rios_km"”), 

B5 = ee.Image("users/lojanjulia/comparacionmodelosIF/dist_antro_km"”), 

B4 = ee.Image("users/incendiosloja/access_vias”), 

ВЗ = ee.Image("users/incendiosloja/access_cen"), 

region =ee.FeatureCollection("users/incendiosloja/loja").geometry(), 
parroquias = ee.FeatureCollection("users/lojanjulia/ParroquiasLoja"), 
listaparroquias = parroquias.aggregate_array(DPA_DESPAR”), 

no_urbano = ee.FeatureCollection("users/incendiosloja/nourbano_provLoja"), 
urbano = ee.FeatureCollection("users/incendiosloja/urbano_cantonLoja"), 
poblados = ee.FeatureCollection("users/incendiosloja/poblados_cantonLoja"); 


( 
( 
(" 
( 


Líneas 228-232: configuración de estilo (tamaño y color de letra) de paneles. 


var с = {}; 
c.directorpanel = ui.Panel([c.info.directorLabel, 
c.info.frreyesLabel],ui.Panel.Layout.Flow('horizontal'), (width: '280px', backgroundColor: 'rgba(255, 255, 
255, 0)'}); 
c.analistapanel = ui.Panel([c.info.analistaLabel, 
c.info.jilojanLabel],ui.Panel.Layout.Flow('horizontal', (width: '280px', backgroundColor: 'rgba(255, 255, 
255, 0)'}); 
c.perfilespanel = ui.Panel([c.info.aboutLabel, c.directorpanel, 
c.analistapanel], ui. Panel.Layout.Flow('vertical'),( 

width: '300px', position: 'middle-left', backgroundColor: 'rgba(255, 255, 255, 0)'}); 
c.infoPanel.add(c.info.TituloMapa).add(c.perfilespanel). add(c.info.websiteL abel); 


Líneas 2-4: llamado a la colección Sentinel-2 


var 525г = ee.ImageCollection(COPERNICUS/S2_SR'); 


var s2Clouds = ee.ImageCollection(COPERNICUS/S2_CLOUD_PROBABILITY5; 


Líneas 955-961: limpieza de la interfaz gráfica e incorporación de los paneles en el mapa. 


ui.root.clear(); 

ui.root.add(c.map); 
c.map.add(c.botones).add(c.infoPanel).add(legena).add(c.controlPanel).add(c.descargaPanel.Title).ad 
d(c.capasPanel).add(button_capasPanel).add(c.descargaPanel);//.add(c.descargaPuntos).add(c.capa 


sPanel); 

c.map.add(c.descargaPanel.chart2); 
c.map.add(c.descargaPanel.chart); 
c.map.centerObject(region).setOptions(ROADMAP”); 
c.dateSlider.setValue(ee.Date(Date.now())); 


Líneas 122-134: configuración y aplicación de parámetros de estilo (tamaño, margen, color) 
en el panel de selección de fecha. 


c.dateSlider = ui. DateSlider(( 


start: "2018-12-13", 
епа: Date.now(), 
value: "2022-08-24", 
period: 5, 


| SELECCIONE UNA FECHA 


se generará el mapa correspondiente 
a la imagen Sentinel-2 anterior más 
1 cercana 
frecuencia: 5 días 


onChange: updateMap, 

style: (width: '180px', position: 'middle-left', 

stretch: 'horizontal', 

padding: 'Opx', + 

margin: '35px Opx Opx Орх', Oct 8, 2022 - Oct 11, 4 
backgroundColor: 'rgba(255, 255, 255, 0)', 08/10/2022 
fontSize: '7px', }, A 
}; 


Líneas 778-788: generación del enlace de descarga para descargar Іа capa de susceptibilidad 
a incendios forestales. 


var URLDescarga = Reclasi_final_MARS.getDownloadURL( ( 
name: 'Susceptibilidad a Incendios Forestales", 
scale: 60, crs:'EPSG:4326', region: region, filePerBana: false, 
format: "GEO_TIFF", }); 
c.descargaPanel. Title.setValue('Capa de Susceptibilidad"); 
c.descargaPanel.Title.setUrl(URLDescarga) ; 


Líneas 610-619: configuración de parámetros de visualización para la capa de límites 
parroquiales y adición al mapa. 


var empty = ee.Image().byte(); 
var outline = empty.paint(( 


featureCollection: parroquias, 
color: 1, width: 1 }); 


var layer4 = ui.Map.Layer( 
outline,(palette: '000000'),'Límites parroquiales', true); 
c.map.layers().set(3, layer4); 


Líneas 357-373: aplicación de filtros (región y fecha) a la colección de Sentinel-2; 
enmascaramiento o “recorte” de nubes y visualización en el mapa. 


function maskClouds (img) {var clouds = ee.Image(img.get('cloud_mask')).select('probability'); var 
isNotCloud = clouds. !t(perminubes); return img.updateMask(isNotCloud);) 
function maskEdges(s2_img) {return s2_img.updateMask( 
s2_img.select('B8A').mask().updateMask(s2_img.select('B9').mask()));) 
var criteria = ee. Filter.and(ee.Filter.bounds(region), ee.Filter.date(START_DATE, END_DATE)); 
var 525 = s2Sr.filter(criteria).map(maskEdges); 
var s2Cloud = s2Clouds.filter(criteria); 
var s2SrWithCloudMask = ee.Join.saveFirst('cloud_mask'.apply([ primary: 525$, secondary: s2Cloud, 
condition: ee.Filter.equals((leftField: 'system:index', rightField: 'system:index'))); 
var s2 = ee.ImageCollection(s2S).median().clip(region); 
var s2CloudMasked = 
ee.ImageCollection(s2SrWithCloudMask).map(maskClouds).median().clip(region); 
var mask_S2 = s2CloudMasked.select('B3'); 
var unmask_S2 = s2.select('B3”); 
var rgbVis = (min: 0, max: 3000, bands: [В4', 'B3', 'B2"]); 
var layer1 = ui.Map.Layer( 
s2, rgbVis, 'S2 without cloudmasking', false); 

c.map.layers().set(0, layer1); 

var layer2 = ui.Map.Layer( 
s2CloudMasked, rgbVis, 'S2 SR masked at' + MAX_CLOUD_PROBABILITY + '%', false); 
c.map.layers().set(1, layer2); 


+ 
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Líneas 412-419: configuración de parámetros de visualización рага Іа capa de registros 
históricos de incendios forestales del SNGRE y adición al mapa. 


var registrosIF = ee.FeatureCollection("users/incendiosloja/Eventos_2010_2021_inf"); 
var style = registrosIF.style({ 
color:'black', width: 1, pointSize: 4, pointShape: 'triangle', fillColor: 'yellow" 


}); 


var layer5 = ui.Map.Layer( 
style,null Registro históricos SNGRE', false); 
c.map.layers().set(4, layer5); 


Líneas 513-545: generación de los índices espectrales de vegetación; transformación de las 
variables a utilizar como “bandas”; unión de las bandas en una sola imagen; aplicación de las 
ecuaciones diferenciales para cada coeficiente en el modelo MARS (para el caso de Loja hay 14 
coeficientes que dividen los datos para su aplicación). 


Hinge function o función bisagra para el Coeficiente 1: 
Coef, = (max (0, (B1 — (-0.0139))) * -127) ‚ siendo B1 = índice МОМ! 


function MARS_Loja (){ var МОМ! = s2CloudMasked.normalizedDifference (['B8A', 'В11']); var B1 = NDMI; 
var NDVI = s2CloudMasked.normalizedDiffterence (['В8', 'B4']); var B2 = NDVI; var bandas = 
B1.addBands(B2).addBands(B3).addBands(B4).addBands(B5).addBands(B6).addBands(B7).addBands(B8).rename('B1''B2','B3''B4','B 
5','B6','B7''B8'); 
var Coef_1 = bandas.expression('float(max(0,(B1-(-0.0139)))*-127)', ('B1': bandas.select('B1'))); var Coef_2 = 
bandas.expression('float(max(((-0.0139)-B1),0)* 82)', ('B1': bandas.select('B1'))); 
var Coef_3 = bandas.expression('tloat(max(0,(B7-1.975))*-3.81)', {'В7': bandas.select('B7'))); var Coef_4 = 
bandas.expression('float(max((1.975-B7),0)*-4.99), {'В7': bandas.select('B7'))); 
var Coef_5 = bandas.expression('float(max(0,(B2-0.3865))*-7.38)', {'В2': bandas.select('B2"))); var Coef_6 = 
bandas.expression('float(max(0,(B7-2.904))*8.02)', {'В7': bandas.select('B7"))); 
var Coef_7 = bandas.expression('float(max(0,(B5-0.303))*0.585)', ('B5': bandas.select('B5'))); var Coef_8 = 
bandas.expression('float(max((0.303-B5),0)*3.54)', ('B5': bandas.select('B5'))); 
var Coef_9 = bandas.expression('float(max((120-B4),0)*0.0485)', {'В4': bandas.select('B4'))); var Coef_10 = 
bandas.expression('float(max(0,(B3-(116)))*-0.653)', {'B3': bandas.select('B3'))); 
var Coef_11 = bandas.expression('float(max(0,(116-B3))*0.619)', ('B3':bandas.select('B3'))); var Coef_12 = 
bandas.expression('float(max(0, (0.1662-B6))*-6.95)', ('B6': bandas.select('B6'))); 
var Coef_13 = bandas.expression('float(max(0,(B3-11))*0.648)', ('B3': bandas.select('B3"))); var Coef_14 = 
bandas.expression('float(max(0,(B1-(-0.0751))*121))', {'В1': bandas.select('B1')); 
var CO = -76.10; var C1= ee.Image(1).where(B1.gt(-0.0139),Coef_1).where(B1.lte(-0.0139),0); var C2= ee.Image(1).where(B1.!te(- 
0.0139),Coet_2).where(B1.gt(-0.0139),0); 
var СЗ = ee.Image(1).where(B7.gt(1.975),Coef_3).where(B7.lte(1.975),0); var C4 = 
ee.Image(1).where(B7.!lte(1.975),Coef_4).where(B7.gt(1.975),0); 
var C5= ee.Image(1).where(B2.gt(0.3865), Coef_5).where(B2.lte(0.3865), 0); var C6 = 
ee.Image(1).where(B7.gt(2.904),Coef_6).where(B7.lte(2.904),0); 
var С7 = ee.Image(1).where(B5.gt(0.303),Coef_7).where(B5.lte(0.303),0); var C8 = 
ee.Image(1).where(B5.!lte(0.303),Coef_8).where(B5.gt(0.303),0); 
var C9 = ee.Image(1).where(B4.lte(120),Coef_9).where(B4.gt(120),0); var C10 = 
ee.Image(1).where(B3.gt(116),Coef_10).where(B3.!te(116),0); 
var C11 = ee.Image(1).where(B3.lte(116),Coef_11).where(B3.gt(116),0); var C12 = 
ee.Image(1).where(B6.!te(0.1662),Coef_12).where(B6.gt(0.1662),0); 
var C13 = ee.Image(1).where(B3.gt(11),Coef_13).where(B3.lte(11),0); var C14 = ee.Image(1).where(B1.gt(- 
0.0751),Coef_14).where(B1.lte(-0.0751),0); var Coef = 
bandas.expression('C0+C1+C2+C3+04+C05+C6+C07+C08+C9+C010+011+C012+013+C14", 

{'Со': CO, 'C1': C1, 'C2": C2, 'СЗ': C3,'C4": C4, 'С5': C5,'C6':C6,'C7':C7, 'C8': Св, 'C9':C9,'C10": C10,C11':011,'C12': C12, 'С13': C13, 
'С14': C14)); 


Se aplica para casos en que el índice NDMI (B1) es mayor a -0.0139; en caso contrario se asigna 0. 
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LANZAMIENTO DEL VISOR COMO UNA EARTH ENGINE APP 


Una vez que el script sea funcional, se procede a revisar el acceso público de las capas o 
variables importadas con el fin de evitar errores y hacer que la experiencia del usuario final 
transcurra sin inconvenientes: 


Share Image: users/lojanjulia/comparacionmodelosIF/Elevacion_div_1000 


lojanjulia gmail.com 


Email or domain Reader + ADD ACCESS 


Select an app v ADD APP ACCESS 


Anyone can read 


Users with access to the asset can view it in the Code Editor via this link: 


https://code.earthengine.google.com/?asset=users/lojanjulia/co f 
mparacionmodelosIF/Elevacion_div_1000 


El lanzamiento del script como un visor de consulta es posible utilizando la opción Apps de la 
esquina superior derecha en nuestro tablero de código. Para una explicación más detallada de 
este proceso se pueden revisar los siguientes enlaces: 


e Publishing an Earth Engine App to Accompany your Research Paper - YouTube 
e Publishing an EE App to Accompany your Research Paper - Documentos de Google 


Publish New App Name and URL 


Name and URL Source Code Publication Details 


App Name and URL 
App Name 


My App 


URL: https://lojanjulia.users.earthengine.app s 


Google Cloud Project €) 
ee-lojanjulia 


CANCEL NEXT 


INSTRUCCIONES DE USO DEL VISOR WEB DE SUSCEPTIBILIDAD A INCENDIOS 


FORESTALES 


https://incendiosloja. users.earthengine.app/view/mars-loja 


El visor web constituye el producto final del proyecto que se ha descrito en este documento 
explicativo. Fue diseñado para ser consultado por todo tipo de usuario: desde los más entrenados 
en Sistemas de Información Geográfica (SIG) hasta miembros de instituciones de respuesta o 
ciudadanos. Por defecto, la página buscará automáticamente la imagen Sentinel 2 más reciente 
disponible, y a partir de ella, mostrará la capa de susceptibilidad a incendios forestales 
acompañada de los límites parroquiales dentro del cantón Loja: 
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En el panel inferior izquierdo, se encuentra el panel de control desde donde se podrá cambiar la 
fecha a consultar de dos maneras: deslizando el cursor de izquierda a derecha (también en 
sentido contrario), o abriendo la opción de calendario. Si se elige una fecha en la que no existe 
imagen satelital disponible, la página buscará la más cercana (recuerde que se actualizan cada 


5 días). 
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А la derecha se encuentra un botón con el texto “Sobre el proyecto”, al hacer clic en él 
se despliega un panel con una breve explicación del proyecto de vinculación, un enlace 
a este documento explicativo y datos de contacto. 
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Al seleccionar el botón “Leyendas” se abrirá un panel con la simbología de las capas de 
susceptibilidad a incendios forestales y de cambio en esta susceptibilidad. 
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Por su parte, el botón “Gráficos de pastel” ubicado еп la parte inferior derecha 
desplegará un resumen visual de la susceptibilidad a incendios forestales, la cobertura 
de nubes y el cambio en el grado de susceptibilidad. 
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Otra función disponible en el visor web es la de acercar el mapa a nuestra ubicación 
actual, siendo mucho más preciso cuando se consulta desde un dispositivo móvil. Para 
regresar al nivel de zoom anterior, bastará con dar clic a “Vista General”. 


25, UTPL s 
58 SmartLand 5 Q Search places Earth Engine Apps 
е Ka [р | “ff kiev я 5 7 | - кии 
ЕГ y Е Kiev Š Dublin Dublin EE 
Е 27 З s zi Suayaquiy ү 
Ss f j ' #2 
А Ру. Е Е 
e j S E Vir 
y 7 А y РА $ Кю ЛӨ та 
г г / o $ y mas E El 
> E 
ES ОД „ЛЕШЕ $ Шү E 
5 гаас = ЕР © © 
Descargar geometrías ё 2 5 > 3 у Ф 
к 3 > ез 3 $ E 
g È Ё 5ш 2 
$ H E 
g 3 g r ] 
z Ibarra $ 
3 ES © 
3 t 2 l Та 
Le é g % + Tulcan " 5 
TEREE = 8 р % 
S gi E WO Лот 
83 FS $. Б А 
3? ASE Guaranda $ A 
с + t y 
SELECCIONE UNA FECHA Ancón C:Ancón | 4e A 
se generará el mapa correspondiente “+ /& $ 
a la imagen БепїпеН2 anterior més q 
CENE ыы soni NL A E 
AN 
Са 
4 19 24 29 Seps 13 18 23 2 Раг EL” A Е 
Н. SA 
|| SENT dl 
# = eira Gamba pe 
2 + Ы 
% Aug 24, 2020 - Aug 29, 2020 RA LS 
24/08/2020 E SÍ Juande salinas * Г 
z Кш! 
¡nue EN > l¿valaivieso— fa Y 
Monteros 7 José Fėlix de Valdivi SDA 
- > ES y у A + © М 
Condo Ф F sel ошо 28.4 "Я 


Para habilitar o deshabilitar las capas de información disponible en el visor es necesario 
desplegar el menú “Layers” de la esquina superior derecha, donde se podrá jugar 
también con transparencias. Finalmente, está disponible la opción de generar 


geometrías (puntos, rectángulos o polígonos) y posteriormente descargarlos en formato 
kml en el menú superior a la izquierda. 
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Mapa Satélite 
Cambios entre el 30/07/2020 y 24/08/2020 
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Focos de calor (7 días) 
Registro históricos SNGRE 
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GLOSARIO DE FUNCIONES DE GOOGLE EARTH ENGINE APLICADAS AL VISOR 
DE SUSCEPTIBILIDAD A INCENDIOS FORESTALES 


ee.Image: un objeto que representa una imagen, generalmente en formato .GeoTIFF. 
Equivalente al uso del término ráster en sistemas de información geográfica tradicionales. 


ee.Image.where(): operador condicional para filtrar o seleccionar pixeles еп ráster (ее.Ітаде). 
Image.where(test, value) 

Ejemplo de utilización dentro del script: 

m.C1= ee.Image(1).where(B6.gt(0.10787),m.Coef_1).where(B6.Ite(0.10787),0); 


En este caso, si la celda tiene un valor mayor a 0.10787 en la banda 6 se reemplaza por el coeficiente 1 
(hipotéticamente valor de 5); si la celda tiene un valor menor o igual a 0.10787 en la banda 6 se reemplaza por 0. 


0.10787 -0.68 0.78 
2.32 1.12 0.012 
-9.85 -5.21 3.25 


ee.Image.byte(): convierte el valor de entrada en un entero de 8 bits sin signo. 
ee.Image.paint(): pinta las geometrías de una colección en una imagen. 
Image.paint(featureCollection,color, width) 
Ejemplo de utilización dentro del script: 
var empty = ee.Image().byte(); 
var outline = empty.paint(f featureCollection: parroquias, 
color: 1, width: 2 }); 

Se utilizó para la visualización de las parroquias en cada cantón sin relleno (límites parroquiales) 

ee.Image.select(): selecciona una banda dentro de un ráster (ee.Image) 
Image.select(var_args) 

Ejemplo de utilización dentro del script: 
var dataset = ee.ImageCollection('FIRMS” filter(ee.Filter.date(fecha2_FC, fecha1_FC)); 
var fires = dataset.select('T21'); 
Con estas líneas de código, seleccionamos solamente una de las tres bandas disponibles en los focos de calor (FIRMS 


Nasa), que es el componente de brightness temperature 


ee.Image.expression(): permite realizar operaciones aritméticas en imágenes ráster, sería 
equivalente a la calculadora ráster en QGIS. 


Image.expression(expression, map) 
Ejemplo de utilización dentro del script: 
m.BSI = s2CloudMasked.expression(float(((Red+SWIR)-(NIR+Blue))/((Red+SWIR)+(NIR+Blue))), 


[Red': s2CloudMasked.select('B11'), 'NIR': s2CloudMasked.select('B8”, 
'SWIR': s2CloudMasked.select('B4”), 'Blue': s2CloudMasked.select('B2”)); 


Cálculo del índice BSI siguiendo la fórmula (Red+SWIR)-(NIR+Blue) / (Red+SWIR)+(NIR+Blue) 
ee.Image.addBands(): integra bandas a una imagen. 


Image.addBands(scrImg, names, overwrite) 


Ejemplo de utilización dentro del script: 


m.bandas = 
B1.addBands(B2).addBands(B3).addBands(B4).addBands(B5).addBands(B6).addBands(B7).addBands(B8).ren 
ame(B1'''B2",'B3','B4"/'B5',B6','B7',B8' ); 


Aquí se integran las ocho variables que se necesitan para el modelo de Quilanga en una sola imagen. 


ee.updateMask(): actualiza la máscara de una imagen en todas las posiciones donde la máscara 
existente no es cero. Equivalente a la herramienta “Cortar ráster por capa de máscara” en QGIS. 


Image.updateMask(mask) 
Ejemplo de utilización dentro del script: 


function maskClouds(img) {var clouds = ee.Image(img.get('cloud_mask').select('probability'); var isNotCloud = 
clouds.It(perminubes); return img.updateMask(isNotClouq);) 


function maskEdges(s2_img) {return s2_img.updateMask( 
s2_img.select('B8A”).mask().updateMask(s2_img.select('B9').mask()));) 


Función para recortar o enmascarar las nubes de las imágenes Sentinel-2 


ee.Image.reduceNeighborhood(): aplica “reductores” a las celdas para suavizar o generalizar 
la imagen ráster. Herramienta equivalente a “Simple Filter” de SAGA GIS. 


ee.Image.tolnt(): convierte el valor de entrada en un entero de 32 bits con signo. 
ee.Image.pixelArea(): permite calcular el área o superficie de cada celda. 
ee.Image.add(): operador matemático de suma (+). 

ee.Image.substract(): operador matemático de resta (-). 

ee.Image.divide(): operador matemático de división (/). 

ee.Image.multiply(): operador matemático de multiplicación (*). 


ee.Image.getDownloadURL() / ee.FeatureCollection.getDownloadURL() : genera un enlace 
para descargar el archivo (ráster o vectorial). El tamaño máximo permitido es 32 MB. 


Image.getDownloadURL (params, callback) 
Ejemplo de utilización dentro del script: 


var URLDescarga = no_nulos.getDownloadURL( { 
name: 'Susceptibilidad a Incendios Forestales, 
scale: 35, 
crsEPSG:32717', 
region: region, 
filePerBand: false, 
format: "GEO_TIFF", 


}) 


En este caso se solicita un enlace para descargar la сара de susceptibilidad a incendios forestales, a una escala de 35 
metros, con el SRC ‘EPSG: 32717" y en formato GeoTiff. 


ee.ImageCollection: un objeto que almacena múltiples imágenes ráster (ee.Image) 
ee.ImageCollection(args) 

Ejemplo de utilización dentro del script: 

var s2Sr = ee.ImageCollection(COPERNICUS/S2_SR”); 


Aquí se hace un llamado a toda la colección de imágenes Sentinel-2 que estén disponibles en Google Earth Engine 
(+5000 imágenes) 


ee.ImageCollection.filterBounds(): permite seleccionar o recortar colecciones de imágenes 
para un límite dado. 


ee.ImageCollection.filterDate(): permite filtrar colecciones de imágenes para un rango de 
fechas. 


ImageCollection.fitterBounds(geometry) 
ImageCollection.filterDate (start,end) 
Ejemplo de utilización dentro del script: 
var criteria = ee.Filter.and(ee.Filter.bounds(region), ee.Filter.date([START_DATE, END_DATE)); 


En este caso se filtra la colección Sentinel 2 para la fecha seleccionada en el panel de control y se recorta para el límite 
del cantón. 


ee.Image.clip(): permite recortar imágenes ráster (ee.Image) siguiendo la geometría de archivos 
vectoriales (ee.Feature / ee.Geometry) 


ee.Image.median(): calcula la mediana de una colección de imágenes. Es posible obtener otros 
estadísticos: .mean(), .min(), .max() 


Image.clip(geometry) 
ImageCollection.median() 
Ejemplo de utilización dentro del script: 
var s2 = ee.ImageCollection(s25).median().clip(region); 
En este caso se calcula la mediana de las imágenes Sentinel 2 para la fecha seleccionada y se recorta para el límite. 


ee.FeatureCollection: un objeto que almacena información de manera similar al formato 
vectorial en QGIS 


ee.Number.parse(): convierte una cadena de texto en número 
ee.Date.fromYMD(): construye una fecha a partir de información base (año, mes, día) 
ee.Date.advance(): crea una nueva fecha según los parámetros que se le asignen 
ee.Date.fromYMD(year,month,day,timeZone) 
Date.advance (delta,unit,timeZone) 
Ejemplo de utilización dentro del script: 


var START_DATE = ee.Date.fromYMD(2020, 08, 24); 
var END_DATE = START_DATE. advance ( 1, 'day”); 


En este caso la fecha inicial sería 24-08-2020, y la fecha final 25-08-2020 (1 día adelante). 
ui.Map.addLayer: añade objetos (ее.Ітаде o ee.FeatureCollection) al mapa como una capa. 
Map.addLayer (eeObject, visParams, name, shown, opacity) 
.layers().set(): permite “enlistar" u ordenar las capas a mostrar en el mapa 
Ejemplo de utilización dentro del script: 
var Іауегі = ui.Map.Layer( 
s2, rgbVis, 'S2 without cloudmasking', false); 
c.map.layers().set(0, layer1); 
var layer2 = ui.Map.Layer( 


s2CloudMasked, rgbVis, 'S2 SR masked at ' + MAX_CLOUD_PROBABILITY + '%', false); 
c.map.layers().set(1, layer2); 


Map.centerObject(): permite centrar el mapa a un objeto dado y a un determinado nivel de zoom. 


Map.setOptions(): permite configurar el mapa base (satélite, terreno, calles, etc.) 

ui.Panel: crea un objeto “widget” o panel que puede contener otros objetos o paneles. 

ui.Label: crea etiquetas de texto que pueden añadirse a los paneles. 

ui.Buttom(): crea un botón para desencadenar una acción. 

ui.Label.setValue(): configurar un valor por defecto para las etiquetas de texto (ui.Label) 
Label.setValue(value) 


Ejemplo de utilización dentro del script: 


c.nubesPanel. Title.setValue('Porcentaje de nubes еп la imagen 527); 
c.nubesPanel.subTitle.setValue('Fecha de la imagen S2”); 
c.nubesPanel.subTitle2.setValue('Permisividad de nubes”); 
c.descargaPanel.Title.setValue('Link de descarga'); 


ui.Select.getValue(): captura el valor elegido en el panel de control para utilizarlo 
posteriormente. 


ui.Button.onClick(): hace una devolución de llamada (repita una acción) cuando recibe un click. 


ui.Select.onChange(): hace una devolución de llamada (repita una acción) cuando cambia la 
selección. 


